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Q '. ABSTRACT 

' Context. The definition of the Galactic coordinate system was announced by the lAU Sub-Commission 33b on behalf of the lAU 

QQ , in 1958. An unrigorous transformation was adopted by the Hipparcos group to transform the Galactic coordinate system from the 

■ FK4-based B 1950.0 system to the FK5-based J2000.0 system or to the International Celestial Reference System (ICRS). For more 
than 50 years the definition of the Galactic coordinate system has remained unchanged from this IAU1958 version. On the basis of 

I I ' deep and all-sky catalogs, the position of the Galactic plane can be revised and updated definitions of the Galactic coordinate systems 

, can be proposed. 

OAims. We re-determine the position of the Galactic plane based on modern large catalogs, such as the Two Micron All-Sky Survey 
' (2MASS) and the SPECFIND v2.0. This paper also aims to propose a possible definition of the optimal Galactic coordinate system 

(~| I by adopting the ICRS position of the Sgr A* at the Galactic center 

Methods. The near-infrared 2MASS point source catalog and the SPECFIND v2.0 catalog of radio continuum spectra are used to 
I . determine the mean position of the Galactic plane on the celestial sphere. By fitting the data to an ideal Galactic equator, the parameters 

O ' defining the Galactic coordinate system are obtained. 

\ Results. We find that the obliquity of the Galactic equator on the ICRS principal plane is about 0.4° (2MASS) and 0.6° (SPECFIND 

C/3 , v2.0) larger than the J2000.0 value, which is widely used in coordinate transformations between the equatorial {a, S) and the Galactic 

{£,b). Depending on the adopted parameters, data, and methods, the largest difference between the resulting Galactic coordinate 
systems is several arcminutes. We derive revised transformation matrices and parameters describing the orientation of the Galactic 
^ ' coordinate systems in the ICRS at the 1 milli-arcsecond level to match the precision of modern observations. 

, Conclusions. For practical applications, we propose that a revised definition of the Galactic coordinate system should be required in 

. the near future. 

■ Key words, astrometry - Catalogs - Galaxy: general - reference system 

(N ■ 
\q ; 

' 1. Introduction ordinates. Here we list the transformation matrix and the pa- 

rameters for the sake of comparison. The transformation matrix 

^ The Galactic coordmate system (hereafter GalCS) is important (3 ^ 3) ^^^^ equatorial to the Galactic coordinate system can 

^-H . for studies of the Galactic structure, kinematics, and dynamics, written as 
sjjl ' In the GalCS, the principal plane (or the Galactic equator) co- 

. !-i ' incides with the plane of the Galaxy and the Galactic longitude = ^3 (90° -6) -Hi (90° - ^p) ■ (90° + a^) , (1) 

^ of 0° is defined as the direction pointing to the Galactic center where (aP, 6^) is the equatorial coordinates of the north Galactic 

; (GC). According to the above criteria, the lAU Sub-Commission pole (NGP) and is the position angle of the Galactic center 

33b summarized previous observations and defined the GalCS (GC) at the NGP with respect to the equatorial pole. These three 

in the framework of the B1950.0 FK4-based reference system parameters were de fined at the epoch of B1950.0 by the lAU in 

in 1958 (Blaauw et al. 1960). The standard GalCS subsequendy 1958 (Blaauw et al. |T960] l. Their numerical values, referred to 

adopted was that for which the pole was based primarily on the the FK5-based J2000.0 reference system, were derived as 

distribution of neutral hydrogen in the inner parts of the Galactic _ . 

disk. Murray (11989) discussed the transformation of the GalCS "^12000.0 - z,u 

from the FK4 to the FK5 system and derived the orientation of ^noooo ~ +27°07'41".704, 

the GalCS at the epoch of J2000.0, and his method was adopted n 1 99° on oi c;-? 

u u- . , A A. f (^12000.0 = 12/ .9319185/, 

by the Hipparcos team as a standard transformation. .op o 

Liu et al. dMB discussed the history of the GalCS, found 'n = 90 - 6^^^^ ^ = 62 .8717, (2) 

that a J2000.0 GalCS which was transformed from the original where /n is the inclination of the Galactic plane. 

IAU1958 definition using unrigorous methods does not match Starting from the J2000.0 system, Liu et al. (1201 II) derived 

the physical features of the GalCS, and pointed out that there the GalCS in the ICRS by using the bias matrix S, and the au- 

may then be some confusions, when applying the Galactic co- thors also tentatively proposed a new definition of the GalCS 
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using the position of the compact radio source Sagittarius A* 
(Sgr A*). They only used a single observation of the GC (ICRS 
position of Sgr A*) to construct a new coordinate system, fail- 
ing to using other data to ensure that this revised definition is the 
optimal choice. Large astrometric and photometric catalogs such 
as the Two-Micron All Sky Survey (2MASS) and other catalogs 
at radio wavelength such as SPECFIND v2.0 have provided us 
with deeper and more homogeneous observational data to esti- 
mate more reliably the position of the Galactic plane. To con- 
struct an optimal GalCS, they also ensured that there was an ac- 
curate match with the real Galactic plane on the celestial sphere, 
namely, that the principal plane {x - y plane) corresponds to the 
median plane of the Milky Way and its x-axis pointing to the 
GC. 

It is known that the Milky Way is a disk-like Galaxy in 
which our Sun is located. The scale height of the thick disk 
is about 1350 pc (Gilmore & Reid 1983) or 1048 pc (Veltz et 
al. 12008 ). Since the vertical displacement of the Sun above the 
Galactic disk (~15 pc) (Zhu 1200^ is much smaller than the 
Galactocentric distance of the Sun (~8.0 kpc) (Reid [1993]) and 
the scale height of the Galaxy, we can assume that the Sun is 
in the middle of the Galactic disk so that the Milky Way is an 
"edge-on" galaxy seen from the Sun, which leads to equal num- 
bers of member stars above and below the mean Galactic plane. 
A deep and homogeneous all-sky survey should reveal the edge- 
on structure of the Milky Way and the stars should be evenly 
distributed above and below a median great circle on the ce- 
lestial sphere, which is the optimal Galactic plane we intend to 
find. In this paper, we select and analyze data from the latest 
near-infrared and radio catalogs to calculate the orientation of 
the optimal GalCS in the ICRS. 

In Sect. 2, we describe the properties of the data that we use 
to fit the mean position of Galactic plane. Sect. 3 presents the re- 
sults for the Galactic plane and the GC positions, and we discuss 
our revised definitions of the GalCS in Sect.4. Finally, in Sect. 
5 we discuss and summarize our main results and draw some 
conclusions. 



2. Data 

We attempt to estimate the position of the Galactic plane using 
data describing the all-sky distribution of stars. Our goal is to 
calculate three parameters (a^, (5^, ff) that describe the orientation 
of the GalCS in the ICRS. A deep all-sky catalog with stars com- 
plete to a certain magnitude would reproduce a smooth Galactic 
belt on the celestial sphere and a highest star density in the di- 
rection of the GC if we did not consider interstellar extinction 
and the vertical distance of the Sun above the Galactic disk. 

A catalog has to meet several criteria if it can be used to 
determine reliably the position of the Galactic plane. First, it 
must cover the area of the whole disk at low Galactic latitudes. 
Secondly, it must be aff'ected by little or no interstellar extinction 
so that the two dimensional (on the celestial sphere) distribution 
of the stars is uniform. Infrared or radio catalogs usually meet 
this criteria but optical catalogs are strongly aff'ected especially 
at low Galactic latitudes. Thirdly, it should contain a sufficient 
number of stars in the framework of the ICRS. The positional 
precision of the individual stars are not crucial because we rely 
on a statistical least squares fit method (see the next section) 
to obtain the most reliable estimate of the mean position of the 
Galactic plane. Proper motions of stars are not needed either, 
because the newly defined GalCS is designed to be epoch inde- 
pendent and fixed with respect to the ICRS. 
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Fig. 1. The sky distribution of the SPECFIND v2.0 catalog in the 
J2000.0 GalCS. The dashed lines are horizontal at b ^ -5° and 
b = +5° 



We choose the 2MASS point source catalog (Cutri et al. 
'2003, Skrutskie et al. '2006) at near infrared wavelengths and 
the SPECFIND v2.0 (VoUmer et al. 2010) catalog in the radio 
band which satisfy all of the above criteria, in our analysis here. 

2.1. 2MASS 

The 2MASS catalog was constructed from uniform observations 
of the entire sky (covering 99.998% of the celestial sphere) in 
three near-infrared bands: J (1.25/im), H (1.65;um), and Kg 
(2.16;um). The limiting flux in each band is about 1 mJy. The 
2MASS contains positions and photometric parameters of 470 
992 970 sources, but no proper motions. The positions of stars 
are calibrated to the Tycho-2 catalog and the astrometric accu- 
racy is 70-80 mas in the magnitude range of 9 < /Tj < 14 mag. 
Figure 13 of Skrutskie et al. (120061) shows a beautiful pictures 
of our Galaxy across the Sky. The 2MASS scans have a typi- 
cal signal-to-noise ratio of ten and is > 99% complete above 
15.8, 15.1, and 14.3 magnitude for J, H, and K,, respectively. 
We choose a complete set of data to fit the Galactic equator by 
limiting the data we consider to those for which 1 1 < 7 < 15.8. 

2.2. SPECFIND v2.0 

The new release of the SPECFIND radio cross-identification cat- 
alog (Vollmer et al. 2010) contains 3.76 million cross-identified 
radio objects, which has about 60% more sources than its 
first release in 2005. The positions of the radio objects in the 
SPECFIND v2.0 catalog were derived from 97 catalogs. Figure 
1 depicts the projected distribution of objects in the J2000.0 
Galactic coordinates system from which we can see clearly the 
Galactic plane, even though the whole catalog is not very homo- 
geneous. In this sense, only objects at low Galactic latitudes are 
selected. 



3. The position of the Galactic plane in the ICRS 

3. 1 . Fitting fhe position of tlie Galactic plane on the celestial 
sphere 

The orientation of the Galactic equator can be defined by 
the equatorial coordinates of the NGP or equivalently 

(ffN, 'n), where otn is the right ascension of the ascending node 
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Celestial equator 

Fig. 2. The pattern of the Galactic plane and the celestial equator: 
cr(a, 6) is a star on the Galactic plane. is the ascending node 
of the Galactic plane with respect to the celestial equator. 



of the Galactic equator on the celestial equator and /n is the in- 
clination of the Galactic plane. They are related by 

aN - = 90°, /n + = 90°. 

Figure 2 shows the relationship between the Galactic plane and 
the equator. The equatorial coordinates (a, 6) of a star cr on the 
Galactic plane satisfies the equation 



tan 6 - sin (a - cn) ■ tan i'n, 



(3) 



in which q-n and /n are two parameters to be fitted from the cat- 
alog data. 

In this paper, we ignore the weak interstellar extinction at 
the infrared and radio bands (which means that the survey ob- 
servation is homogeneous) so that the Galaxy has a symmetrical 
structure as seen from the Sun. The optimal Galactic equator 
should be a great circle on the celestial sphere that is the clos- 
est fit to the density distribution pattern of the above catalogs: 
the position of the Galactic plane is where the highest star den- 
sity appears in the direction of the Galactic longitudes, because 
stars should be concentrated on the real Galactic plane due to the 
gravitational attraction and rotation of the Milky Way. In other 
words, the great circle chosen should have an equal number of 
stars on both sides. 

For the 2MASS point source catalog, we select stars within 
the magnitude range 11 < J < 15.8. The upper limit to the 
magnitude range is the 99% completeness magnitude in the J 
band as mentioned in Sect. 2.1, and the lower limit rejects bright 
stars close to the Sun. Owing to the local structure of the of the 
Galaxy, such as the Gould belt (Westin 1985 ) and the Galactic 
warp (Miyamoto et al. 1988), which may affect our results, we 
exclude nearby stars by applying a bright upper limit to our se- 
lected sample. The complete sample comprises about 56.7 per- 
centage of stars in the whole 2MASS catalog. The density map 
(360 bins in f x 1 80 bins in /? ) of the selected stars in the J2000.0 
Galactic coordinate system is shown in Fig. 3, which depicts a 
beautiful edge-on Galactic plane on the sky. The star density at 
high Galactic latitudes is much lower than closer to the Galactic 
plane and the influence of halo stars is not high. 

Instead of fitting Eq. (3) to individual star coordinates (a, i5), 
we choose to use the equatorial positions of the geometrical cen- 
ters of each 1 "-Galactic-longitude bin (360 in total), in order to 
derive a more accurate estimate of ckn and /n. We analyze the 
2MASS data in this way for the following reasons. The stellar 
density close to the GC is extremely high (which means that 
the weight of the data in the GC region is much larger than that 
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Fig. 3. The star density distribution of the 2MASS point source 
catalog in the magnitude range 11 < J < 15.8. The direction of 
GC (^ = 0, ^7 = 0) is plotted at the center 

in any other area), such that the small-scale structures close to 
the GC will strongly affect the resulting position of the Galactic 
plane. Our processing method assigns to data at all Galactic 
longitudes the same weight to avoid any biasing. The densities 
within the LMC and SMC areas are excessively high in the cur- 
rent study and are replaced with suiTounding densities in subse- 
quent calculations. 

We first calculate the central Galactic coordinates of the 
360 columns of the 360x180 density-distribution matrix and 
then transform (4, bk), k - 1 • ■ • 360 into an equatorial co- 
ordinate system to derive the equatorial coordinates {ak,Sk), 
k - 1 • ■ • 360, by applying the previous J2000.0 transformation 
matrix in Eq. (1). In the next step, these equatorial coordinates 
are fitted with Eq. (3) to obtain otn and for the first time. With 
the additional position angle 0j2ooo.o in Eq. (2), we can write the 
new transformation matrix N, according to which the updated 
Galactic coordinates of the stars and the updated density matrix 
of our samples can be obtained. Repeating the above processes, 
we derive converged parameters aN and /n, which have the nu- 
merical values 

OTN = 282°.987 + 0.032, /n = 63°.280 + 0.019. 

These values of {a^, /n) are equivalent to the lAU definition pa- 
rameters (q'P,5P) and we can quote them to infinite accuracy 
(here and hereafter to 1 milli-arcsecond), if they are adopted as 
a definition 



'2MASS 



12"51"56'.7726, 



^2MASS = +26°43'11".096. 



(4) 



It should be understood that this is a numerical convention and 
does not necessarily imply that any given parameters and ele- 
ments of the matrices (in the following equations) are known to 
the quoted accuracy. The standard 1 - cr errors in above results 
reflect only the uncertainty in the data fitting and have no deep 
physical meanings. 

Figure 4 illustrates the positions of the 360 bin centers, in 
which the gray curve shows the original value, whereas the black 
one is plotted after and /n converge. The raw and reduced 
curve are almost evenly situated above and below b - Q. They 
have similar shapes along the great circle b - Q and correspond 
to the old and new transformation matrix N. The difference be- 
tween the two curves reflects the small-angle orientation bias 
between the original and the final GalCSs. 

In the SPECFIND v2.0 catalog (see Fig. 1), objects are not 
evenly distributed on the sky: the southern hemisphere was ob- 
served less than the northern. The global structure of this distri- 
bution has a strong effect on our results, even though the Galactic 
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Galactic longitude (deg) 

Fig. 4. The curve of the center in 1° longitude bins. The gray 
curve is the original curve transformed by J2000.0 rotation ma- 
trix and the black one the reduced curve transformed by new 
matrix that is derived from 2MASS fitting. 

structure can be clearly seen in Fig. 1 . We only retain the data 
at low Galactic latitudes \b\ < 5°, which avoids possible biasing. 
We divide the belt (0° < Z < 360°, -5° < ^7 < 5°) into 360 bins, 
each one with l°(in £) x 10°(in b) coverage. The mean positions 
of objects in each bin (^k, ^k), ^ = 360 are then calculated 

and transformed into the J2000.0 equatorial coordinates (Sk, 5k), 
k - l,---,360. Similarly as we process the 2MASS catalog, 
these 360 equatorial coordinates are used to to fit au and Zn in 
Eq. (3) using a least squares method. After we obtain the result- 
ing Q-N and /n, the equatorial to Galactic transformation matrix 
N can be derived by applying 6'j2ooo.o- The equatorial coordinates 
of the raw data can then be transformed into the Galactic coor- 
dinates using the new matrix. In the next step, we recalculate the 
updated coordinates of the mean positions and fit Eq. (3) 

to get another pair of and /n values. These procedures are 
repeated until the two objective parameters converge to the final 
values. The resulting parameters are 

an = 282°.699 + 0.031, /n = 63°.513 + 0.020, 

or the coordinates of the NGP are 



«SPECFIND - 12''50-47^6503 
-^SPECFIND = +26°29'12".412. 



(5) 



We note that the position angle 0j2O()O.o (which indicates the 
position of the GC on the Galactic plane) is unaltered in itera- 
tive steps because we are only concerned with the position of 
the Galactic plane in both the 2MASS and SPECFIND v2.0 data 
sets. From Eqs. (4) and (5), we find that the fitted values of in- 
clination /n from the 2MASS and the SPECFIND v2.0 are about 
0.41° and 0.64° respectively, larger than the J2000.0 reference 
values in Eq. (2). The ascending node fitted using the 2MASS 
data is about 7 arcseconds eastward of the J2000.0 node, but the 
corresponding result using the SPECFIND v2.0 data is about 10 
arcseconds westward of the J2000.0 node. We emphasize here 
that our results depend on the choice of data (range of magni- 
tude, type of stars, bandwidths, etc.), although both 2MASS and 
SPECFIND v2.0 are appropriate candidate catalogs to fit and 
iN because of their advantages described above. The generated 
results in Eqs. (4) and (5) present, to a certain extent, the po- 
sition of the Galactic plane in the near-infrared and radio con- 
tinuum data sets based on modern observations. We confirm that 
complete and deeper long-wavelength surveys will provide more 
reliable estimates. Although there are many extremely large op- 
tical catalogs, suc h as U CAC3 (Zacharias et al. ISOTOl l. USNO 
Bl.O (Monet et al. 120031 . and the latest PPMXL catalog (Roeser 



et al. 1201 Oi l (of nearly one billion stars), they cannot be used in 
this study, mainly because of the non-isotropic interstellar ex- 
tinction along different lines of sight, which destroys the homo- 
geneity of the density distribution of stars in these optical cata- 
logs. 

3.2. The position of tlie Galactic center in the ICRS 

The lAU working group reviewed the precise position of the 
compact radio source Sagittarius A*, in the B1950.0 reference 
system (Gum & Pawsey 1960) to be 



'B1950.0 



B 1950.0 



-28°57', 



which was afterward used to define the direction of the x-axis 
(GC) of the IAU1958 GalCS. According to the rigorous FK4- 
FK5 transformation proposed by Standish ( 1982j and Aoki et al. 
(il983j . the position of the Sgr A* at J2000.0 is 



SgrA* _ , -h.jrm.^s 7 cSgiA. 

"J2000.0 ~ ^ ' ^-^ ^ ' • ' ' "jaooo.o 



-28°.58'09" 



(6) 



while the x axis direction of the J2000.0 GalCS recommended 
by Hipparcos team (ESA |1997] l is 

<ooo.o = 17H5™37M991, ^"jooo.o = -28°56'10".221. (7) 

Using Very Long Baseline Array (VLBA) measurements, Reid 
& Brunthaler (2004) derived the absolute position (referred to 
the ICRS) of the Sgr A* at the dynamical center of the Galaxy 

o'Reid* = 17''45"'40'.0400, ^dd' = -29°00'28".138. (8) 

Since the Milky Way is a gravitational system centered at the 
GC, the stellar density should be highest in that region. We find 
from the 2MASS data set that the position of the peak density is 



^^2MASS = 17*>46'"6^5328, 5° 



2MASS 



-28°54'16".164. 



(9) 



In the following section, we adopt the Sgr A* position of Reid 
& Brunthaler (2004) derived using direct radio observations as 
the direction of the GC because it is independent of complicated 
transformations. 



4. Constructing the Galactic coordinate system 

We denote the orthogonal GalCS using three principal axes 
[x,y,z]. The orientation of the IAU1958 and Hipparcos 
(J2000.0) GalCSs are characterized by the coordinates of the 
NGP (z-axis) and the position angle of the GC (jc-axis). After 
this set of parameters has been determined, we can construct 
the GalCS in the corresponding equatorial reference system us- 
ing the 3x3 matrix N. In Sect. 3, we have derived possible 
the positions of the NGP (Galactic plane) (aP,5P) and the GC 
(ff", (5°). However, these coordinates or angles cannot be used to 
construct the GalCS directly, because their directions in (q-p, i5P) 
and (ff", (5°) were determined independently, either by catalog 
fitting or by observations, and they are not perpendicular to each 
other. In other words, a given GC may not lie on a given Galactic 
plane. We construct the GalCS in two diff'erent ways: (1) we as- 
sume that the z-axis coincides with the given direction of the 
NGP (ffP, (JP) then try to find the x direction using the given GC 
coordinates (0-°, <5'^); (2) we assume that the x-axis coincides with 
the given direction (o-", 5'^) of the GC then try to infer the orien- 
tation of the Galactic plane from catalog fitting. 
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4.1. Method 1: z-axis fixed condition 

We first assume that the z-axis of the GalCS points in the direc- 
tion of the NGP (q"" = q-p, 6^ = 6^) that is fitted from either 
the 2MASS or the SPECFIND v2.0 catalog. Since the given GC 
does not lie on the Galactic plane (x-y plane), we try to find the 
nearest point (a", 6'') from the given GC (a", 5'^) on the fixed 
Galactic plane (q;P, S^). Evidently, {a-\ 6^) and (a", 5") share a 
common position angle 6 with respect to the north celestial pole 
as seen from the NGP, given by 



sm 



tan 6* 



cos 5P ■ tan (5" - sin 6^ ■ cos (a^ - aP) 
Substituting in both Eqs. (4) and (8), we find that 
OiMASS = 123°.044663150. 



(10) 



(11) 



On the basis of these values of a'2MASS' ^2mass' ^2mass, we 
derive our new equatorial-Galactic transformation matrix based 
on the 2MASS data set and the Sgr A* position (printed to 10 
decimals, 1 milli-arcsecond accuracy) 

A/2MASS = 

-0.0505347007 -0.8719028362 -0.4870643574 ^ 
+0.4897972973 -0.4466482209 +0.7487349160 . (12) 
-0.8703705255 -0.2007257110 +0.4496268868, 

If we adopt Eqs. (5) and (8) to compile the GalCS, the posi- 
tion angle of x-axis and the rotation matrix are 

0SPECFIND = 122°.914360525, (13) 
and 

A/SPECFIND = 

-0.0518807421 -0.8722226427 -0.4863497200 ^ 

+0.4846922369 -0.4477920852 +0.7513692061 . (14) 

-0.8731447899 -0.1967483417 +0.4459913295, 

4.2. X-axis fixed condition 

In the above calculations, the z-axis of a constructed GalCS co- 
incides perfectly with the fitted NGP, but the x-axis direction 
does not match the given GC (Sgr A*), but instead offset by sev- 
eral arcminutes. We now assume that the x-axis coincides with 
the given GC (a^ - oq, S" - <5o), and determine the coordi- 
nates of the NGP (z direction) by fitting either the 2MASS or 
SPECFIND v2.0 catalogs with the same strategy as in Sect. 3. 
Since we have already determined two of the three parameters 
that define the GalCS, there remains only one parameter to be 
determined. We choose the position angle rj of the NGP (seen 
from the given GC) as the last unknown parameter. We refer to 
Fig. 5 for the definition of rj. The J2000.0 reference value of rj is 

'7J2000.0 = 58°.599, 

and the equation for fitting is then 

cos (5° • tan<5 = sin<5'' ■ cos^c - o-*') + sin^a - c") ■ tan 77. (15) 

The equatorial coordinate system to GalCS transformation 
matrix using {a^, S^, rf) has the form 



N^Kx {!]) ■ Hi (-5°) • ^3 («°) . 



(16) 



With the matrix N, three standard parameters a^, S^, and 9 of 
the lAU scheme can be determined. In particular, we specify the 




Fig. 5. The position angle 77 of the NGP at the given GC (a°, SP), 
where 9 is the position angle of the GC, cr(a, 6) is a star on the 
Galactic equator, and P is the celestial pole. 



coordinates in Eq. (8) to be the x-axis and fit the 2MASS stellar 
density data iteratively with Eq. (15) to determine rj and N. The 
final results are 



'72MASS =58°.947 + 0.025, 



A/2MASS = 

-0.0546572359 -0.8728439269 -0.4849289286 ^ 
+0.4888603641 -0.4468595864 +0.7492209651 
-0.8706481098 -0.1961121855 +0.4511229097, 



«2MASS = 12''50'"4615444, 



(17) 



(18) 



''2MASS 
^2MASS 



+26°48'56".706, 
122°.912729244. 



(19) 



Similarly, we obtain the orientation of the GalCS based on 
Eq. (8) and the SPECFIND v2.0 catalog 



/7SPECFIND = 59°.275 + 0.024, 



A/SPECFIND - 

-0.0546572359 -0.8728439269 -0.4849289286 ^ 
+0.4838685275 -0.4479748647 +0.7517910405 
-0.8734322153 -0.1935510264 +0.4468267735, 



(20) 



(21) 



' SPECFIND 



^2h49m3gs^73gQ 



^SPECFIND = +26°32'24".989 
0SPECFIND = 122°.823292026. 



(22) 



The resulting values of rj from the both 2MASS and the 
SPECFIND v2.0 differs by a few tenths of a degree when re- 
ferred to 77J2000.0- 

The GalCSs derived by applying the two methods with the 
same observational data (2MASS or SPECFIND) differ in terms 
of the small-angle orientation by several arcminutes, while the 
bias in the GalCSs caused by the use of different observational 
data is about a few tenths of an arcsecond. In Sect. 4 of Liu et al. 
(1201 11 1, the authors proposed a possible definition of the GalCS 
based on the ICRS position of Sgr A* and the transformed 
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GalCS, while the present proposals depend only on the obser- 
vational data but are unrelated to the old GalCS and its transfor- 
mations between various fundamental reference systems, which 
should be avoided when defining a new coordinate system. In 
the context of our aim to establish a revised GalCS, we have 
succeeded in defining an optimal GalCS in which the z-axis co- 
incides with the NGP and the x-axis with the position of Sgr 
A*, which is regarded as the GC. All data have been used were 
measured in the framework of the ICRS. The GalCS should be 
independent of epoch, thus the transformation matrix should be 
applied at the epoch of J2000.0. The results presented in Eqs. (8), 
(11), (13), (19), and (22) provide possible references for future 
revisions of the definition of the GalCS. 
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5. Discussion and conclusion 

On the basis of near-infrared and radio observational data, we 
have described several approaches for calculating the positions 
of the Galactic plane and the GC. From our fitted results, we 
have found that the inclination of the Galactic plane with respect 
to the ICRS equator is either 0.4° (2MASS) or 0.6° (SPECFIND 
v2.0) larger than commonly used values recommended by the 
Hipparcos. This impfies that there remains room for improve- 
ment in the J2000.0 GalCS. 

As mentioned at the beginning of the paper, the main pur- 
pose of establishing a robust GalCS is to ensure that in studies 
of Galactic kinematics, dynamics, and structure, the Galactic co- 
ordinates ({, b) of an object are always reUably and universally 
transformed from its equatorial coordinates {a, S). We note that 
the establishment of the GalCS need not have the same degree of 
precision as either the ecliptic or equatorial reference systems. 
The FK4, FK5, and Hipparcos or the ICRS are reference sys- 
tems based on carefully determined positions of the equator and 
equinox, while the Galactic coordinate system is based on the 
concentration of tens of millions of stars or continuum radio ob- 
jects in a statistical sense. We have defined in the previous sec- 
tions a transformation method (aP, 6^, 6) that is the most appro- 
priate choice for Galactic studies, at least at near-infrared and ra- 
dio wavelengths. This transformation definition should be spec- 
ified in certain fundamental reference systems, the ICRS. Our 
philosophy is similar to the IAU1958 method and the results are 
given for two bands, both based on the latest observations. The 
parameters presented in Sect. 4 provide some possible formal 
definitions of the GalCS that may be applied in the future. 

With the development of astrometry to an unprecedented ac- 
curacy and the acquisition of a large amount of data in recent 
years, we recommend that the Galactic coordinate system be de- 
fined based on observational data rather than a transformation of 
the original definition. The optimal definition of GalCS remains 
an active field of research that should be considered by relevant 
lAU commissions. 
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